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ABSTRACT 

We report on radio observations of five magnetars and two magnetar candidates carried out at 1950 
MHz with the Green Bank Telescope in 2006-2007. The data from these observations were searched 
for periodic emission and bright single pulses. Also, monitoring observations of magnetar 4U 0142+61 
following its 2006 X-ray bursts were obtained. No radio emission was detected was detected for any 
of our targets. The non-detections allow us to place luminosity upper limits of iigso ^ 1.60 mJy kpc^ 
for periodic emission and ^1950, single ~ 7.6 Jy kpc^ for single pulse emission. These are the most 
stringent limits yet for the magnetars observed. The resulting luminosity upper limits together with 
previous results are discussed, as is the importance of further radio observations of radio-loud and 
radio-quiet magnetars. 
Subject headings: pulsars: individual 



1. INTRODUCTION 

Magnetars are a subclass of pulsars powered by the 
decay of their ultra-st rong magnetic fieldSj typic ally 
B ~ io W^ G ^see lWoods fc Thompson! (pOM ). or 
iMereghettil (|2008[ ) for recent reviews). The X-ray emis- 
sion of magnetars exhibits phenomena not seen in X-ray- 
detected rotation-powered pulsars. Magnetars typically 
have a greater X-ray luminosity than can be explained 
by their spin-down luminosity (i.e. Lx > E), unlike 
rotation-powered pulsars. In addition, magnetars are 
very variable objects; they have a wide variety of X-ray 
emission behavior such as short single (~ 100 ms) bursts, 
periods of outburst containing many short bursts, giant 
flares lasting hundreds of seconds and which exhibit pul- 
sations in their fading tails, flux enhancements lasting 
hundreds of days and X-ray pulse profi l e variations (see 
Woods fc Thompson|[2006l: iKaspil [20071: iRea fc Espositol 
201 IL for more details) . Magnetars also have complicated 
timing properties such as timing nois e, torque variations , 
and glitches (e.g . Kaspi ct al. 2000; D an'Osso eralll2003l : 
IDib et al.l[2007ll2008l: |Dib|l2009) . 

According to the McGill SGR/AXP Online Catalog, 
there are currently 16 confirmed magnetars, consisting of 
seven Soft Gamma Repeaters (SGRs) and nine Anoma- 
lous X-ray Pulsars (AXPs), as well as seven magnetar 
candidates (four SGR and three AXP candidates). These 
numbers are small relative to the total number of known 
rotation-powered pulsars, near ly 2000, according to th e 
online ATNF Pulsar CatalofI (jManchester et al.ll2005[) . 

Prior to 2006, there was no firm detection of pulsed 
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radio emission from a magnetar. Since that time, 
three magnetars have been found to have emission at 
radio frequencies: XTE J 1810-197, IE 1547 .0 -5408 
and PSR J1622-4950 (|Camilo et all [20061 I2007bl: 
ILevin et al.|[20To[ ). The radio emission of these three 
sources has properties that are common among them, 
but which are not shared with the bulk of the rotation- 
powered radio pulsar population. Specifically, the three 
radio magnetars show variable spectral indices, vari- 
able pulse profiles, and variable radio luminosities. The 
spectra observed from radio-detected magnetars (even 
when taking into consideration the variability observed 
in their spectral indices) are generally very fiat, or ris- 
ing with observing frequency, in contrast to the majority 
of rotation-powered pulsars, whose spectra steeply de- 
cline with observi ng frequency, with m ean spectral index 
(a) = -1.8 ± 0.2 (jMaron et al.|[2000l ). The origin of the 
spectral flatness in radio emission from magnetars is not 
known. 

The discoveries of radio-loud high-magnetic-field pul- 
sars and magnetars have spurred observations of other 
magnetars in the hope of detecting ra dio emission 
(jBurgav et al.l 120061: [&awford et all I2007D . Limits of 
'S'i4oo ~ 20 /xJy have been placed on three South- 
ern magnetars, IRXS J170849.0-400910, IE 1841-045 
and IE 1048.1—5 937, and one mag n etar candidate, 
AX J1845-0258 (IBurgav et al.l [20Qa ICrawford et all 
[2001 . Also, upper limits on the flux density of sin- 
gle pulses at 1400 MHz from the above-listed four 
sources have been p laced in the range 0.9 - 1.1 Jy 
([Crawford et al.l[2007[) . 

In this paper we describe similar radio observations of 
magnetars in the Northern sky. The goal was to dis- 
cover more magnetars that exhibit pulsed radio emis- 
sion. Increasing the number of radio-loud magnetars 
known could offer insight into the pulsar emission mech- 
anism, the behavior of matter in ultra-strong magnetic 
fields, as well as the possible evolutionary relationship 
between magnetars and the much larger known popula- 
tion of rotation-powered pulsars. 

This paper is organized as follows: ^ describes the 
observations undertaken. The analysis performed on the 
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data collected is described in 331 The results are pre- 
sented in 21 In SjSlwe put the results into context and 
discuss their implications. 

2. OBSERVATIONS 

Observations of five confirmed magnetars and two 
magnetar candidates were carried out using the NRAO 
100-m Green Bank Telescope0 in 2006-2007. The goal 
was to observe the magnetars and magnetar candidates 
visible from the GBT in order to detect radio emission, 
or in the absence of a detection, establish baseline mea- 
surements should a target turn on as a radio source some 
time in the future. The proposal also included Target- 
of- Opportunity (ToO) observations to be triggered if a 
source exhibited an outburst, as determined by on-going 
X-ray monitoring observations. 

Total intensity data were recorded by the GBT's Pul- 
sar Spig ot backend , an auto-correlation spectrometer 
(|Kaplan et al.ll2005l) . Lagged-products were converted to 
spectra off-line. The result is spectra containing 600 MHz 
of usable bandwidth centered at 1950 MHz (S-band) di- 
vided evenly into 768 channels, written out every 81.92 

IIS. 

2.1. Targets Observed 

Here we present relevant details of our seven targets. 
In total, 19 observations of a variety of durations were 
made between 2006 November and 2007 October as sum- 
marized in Table [H For each source we present its best 
distance estimate, together with an estimate of the free 
ele ctrons along the l ine-of -sight toward the source from 
the lCordes fc Laziol ()2002[ ) model. The free electron con- 
tent is parameterized by the dispersion measure (DM). 
Estimates of the model are used to determine the upper 
limit on DM searched for each target. In all cases the 
DM searched was >2 times larger then the maximum 
DM predicted along the line-of-sight. 

AXP IE 18 41—045: This source ha s a rotational pe- 
riod of 11.8 s fjVasisht &: Gotthelj [19971) . The estimated 
distance to this source is kpc (ITian fc Leahvl 

I2008D . which, together with its position, suggests a DM 
~ 800 cm~'^ pc (ICordes fc Laziol [200l . The supernova 
remnan10 Kes 73 is associated with IE 1841-045. The 
pulsar is monitored regular ly by NASA's R ossi X-ray 
Timing Explorei^ {RXTE, (Dib et al.|[2008H Dib 2009). 
The analysis of these data allow for an accurate rota- 
tional ephemeris for the pulsar to be determined for our 
GBT observing epoch; this is shown in Table [21 

IE 2259+586: Th e magnetar has a period of 6.98 s 
([Fahlman fc Gregorvl [T981'). The source is estimated to 
be at a distance of 4.0 ± 0.8 kpc ()Tian et al.ll2010f ). This 
distance and the source position correspond to a line-of- 
sight DM ~ 150 cm-3 pc. IE 2259+586 is associated 
with the supernova remnanlQ CTB 109. IE 2259+586 
is also monitored regularly using RXTE (Dib 2009,), and 
the resulting ephemeris is shown in Table [2l 

* http;//www. gb.nrao.edu/gbt/ 

^ Kes 73 is also referred to as G27. 4+0.0 in Green's 
online Catalog ue of Galactic Super nova R emnants: 
|http: / /www . mr ao. cam .ac.uk/ surveys /snrs/ ( GreenI 120091 ) . 

Information about RXTE monitoring targets and 
observations can be found on the telescope's website: 
littp: / /heasarc. gsfc.nasa.gov/docs/xte/. 

CTB 109 is also known as G109. 1-1.0 in Green's catalogue. 



4U 0142+61: This source has a period of 8.69 s. The 
distance to this source is estim ated to be 3.6 ± 0.4 kpc 
(jDurant &: van Kerkwiikll2006[ ). This distance suggests 
a DM c=L 100 cm~^ pc. The rotational ephemeris for 
4U 0142+61 for our observing epoch is shown in Table 
[2l This ephemeris was also computed using RXTE da ta 
from regular monitoring observations ()Dib et al.ll2007[ ). 

Within a few months before and after the start of 
our GBT observations. X -ray bursts were detected from 
4U 0142+61 with RXTE (jGavriil et al.|[2007l[20TT[ ). The 
six bursts were detected in three separate observations on 
MJDs 53831 (one burst), 53911 (four bursts) and 54138 
(one burst). The last of the six bursts had the largest 
peak X-ray flux. Monitoring observations of 4U 0142+61 
were made on the day of this largest burst, as well as at 
1 day, 1 week, 3 weeks, 2 months and 8 months following 
this burst (see Table [1]). 

SG R 1806—20: The source h as a rotational period of 
7.56 s (jKouveliotou et aLlll998| ) and is estimated to be 
at distance of 8.7;j J kpc (|Bibbv et al.ll2008[ ). which sug- 
gests DM ~ 750 cm"3 pc. SGR 18 06-20 was observe d 
to have a giant flare in late 2004 (jHurlev et all I2005D . 
Following the flare, a fading, extended, impulsed radio 
nebu la was observed around the source (jGaensler et al.l 
120051 ). Two days after the flare, a search for pulsed radio 
emission was performed at 1400 MHz using the Parkes 
Radio telescope. No pulsed e mission was detected down 
to a limit of 51400 ^ 0.2 mJy (jGaensler et al.|[2005l ). 

SGR 1900+14: This S GR has a periodicity of 5.17 
s (jKouveliotou et aLlll999f ). Its distance is estimated to 
be 12 - 15 kpc, which yields an estimated DM ~ 700 
cm~^ pc, using the upper limit of the distance estimate. 
SGR 1900+14 was observ ed to have a giant flare in mid- 
1999 (jHurlev et al.lll999[) . Unpulsed, fading radio emis- 
sion was also observed ar ound SGR 1900+14 following 
the flare (jFrail et al.lll999f ). 

AX J1845— 0258: This magnetar candidate was 
detected in a single 1993 Advanced Satellite for Cos- 
mology and Astrophysics (ASCA) observation wit h 
P = 6.97 s (iTorii et al.l 119981: IGotthelf fc Vasishtl [19981) . 
AX J1845— 0258 is tentatively classified as an AXP can- 
didate because of its long period, soft spectrum, and 
low luminosity. The distance to the pulsar is estimated 
to be 8.5 kpc (jTorii et al.lll998l ). This distance corre- 
sponds to a DM ~ 750 cm-^ pc. AX J1845-0258 is spa- 
tially coincident with the supernova remnant G29. 6+0.1 
(jGaensler et al.l[T999t ). suggesting that the remnant and 
magnetar candidate are associated. N o subsequent de- 
tcction of this pulsar has been made fsee lTam et al.|[2006l 
for example) , hence we have only an approximate period 
estimate. 

GRB 050925: This 7-ray burst (GRB) was detected 
as a short burst (104 ms) with the Swift s atellite's burst 
alert telescope (BAT) on Se pt. 25, 2005 ([Holland et al.l 
[200I ISakamoto et"al1 [20T1I) . GRB 050925 is consid- 
ered an SGR candidate because of its soft X-ray spec- 
trum, which is uncharacteristic of short GRBs, and its 
localization in the Galactic plane {I — 72.32(3)°, h = 
—0.09(3)°). Within the BAT error circle two variable 
X-ray sources are plau sible counterparts of GRB 050925 
([Sakamoto et al.ll2011[ ). GRB 050925 was also foUowed- 
up at 4900 MHz with the Westerbork Synthesis Radio 
Telescope. No sources were detected within the BAT 
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error circle with 5*4900 > 76/xJy, and one constant ra- 
dio source was dete cted just outside the error circle 
(|van der Horsli2005hl lal). This magnetar candidate is still 
unconfirmed, without a period measurement, or an es- 
timate of its magnetic field strength. Nevertheless we 
observed the position of GRB 050925 to search for ra- 
dio bursts or pulses from this possible magnetar. For 
the purpo ses of this work w e use the distance estimate 
from Sak amoto et al.l (|201lD . 8.8 kpc, which assumes 
GRB 050925 is associated with the Hn region W 58. This 
distance and lin e-of-sight corresp ond to a DM estimate 
of 300 cm-3 pc (iCordes &: Laziol[200l ). 

The properties described above are summarized in Ta- 
ble! 

3. DATA ANALYSIS 

Data were searched for both periodic signals and for 
bright individual pulse s using the PR ESTO suite of pul- 
sar search program^ (|Ransomll2001l ). Masks to remove 
narrow-band impulsive and periodic radio frequency in- 
terference (RFI) were created using PRESTO 's rfifind 
program. The masks produced by rf if ind were adjusted 
slightly by hand to mask any additional bad frequency 
channels or time intervals from the data. These masks 
were applied to the data before dedispersion. 

In seven data sets, the mean level had a large step 
during the observation. This was a result of large 
bursts of broadband RFI, which caused offsets in the 
mean level of GET Pulsar Spigot data (S. Ransom, 
private communication). This behav ior has also been 
seen in other Spigot data sets (e.g. iKondratiev et al.l 
120091 at 820 MHz). Our observations that were af- 
fected are: 4U 0142-^61 (MJDs 54056, 54059, 54146), 
IE 1841-045 (MJD 54052), IE 2259-f586 (MJD 54059), 
and SGR 1806-20 (MJDs 54190, 54211). To deal with 
these jumps we did the following. If the jump in the 
data mean occurred near the beginning or end of an ob- 
servation, only the longer portion of the observation was 
analysed. The integration time of the observation was 
adjusted accordingly. Alternatively, if the jump occurred 
near the middle of an observation, the observation was 
split into two portions and each was analysed indepen- 
dently. 

Data were dedispersed according to the following plan. 
For DMs in the range to 1247 cm"'' pc, DM step sizes 
of ADM = 1 cm~^ pc and downsampling by a factor of 
4 were used. In the range DM = 1248 to 2446 cm~^ pc, 
data were downsampled by a factor of 8 and the ADM 
used was 2 cm"'' pc. Finally, for DMs larger than 2448 
cm~^ pc, ADM — 3 cm"'' pc was used and data were 
downsampled by a factor of 16. Maximum DM values 
used varied by target, depending on the p rediction for the 
given line-of-sight by the NE2001 model ( Cordes &: Laziol 
I200I see Table [3] shows the maximum DM values 
searched for each target. Dedispersed time series were 
produced using PRESTO 's prepsubband. 

The resulting time series were inspected for steady 
periodic signals using PRESTO's accelsearch, which 
searches for peaks in Fourier transforms. Sensitivity to 
narrow pulses was increased by summing up to 16 har- 
monics. The search for periodic signals was repeated 
with a red- noise removal technique, which included sub- 

^'^ http;/ /www. cv. nrao.edu/~sransom/presto/ 



traction of a running median from the time series. 

For the three magnetars with known rotational 
ephemerides (see Table [2|) , the data were folded at the 
known ephemeris to search for periodic emission. Even 
when an ephemeris was available, the dedispersed time 
series were still searched using accelsearch for other 
periodic sources in the field-of-view. 

Significant periodicity candidates were folded using 
the DM trial for which the candidate's period had the 
largest signal-to-noise ratio. Folding was performed us- 
ing prepf old. The resulting plots were examined by eye. 

Bright single pulses were searched for using PRESTO's 
single_pulse_search.py, which uses matched filtering 
with top-hat filters having widths ranging from 0.33 to 
50 ms. The data were then downsampled by a factor of 
4 and re-searched to increase sensitivity to single pulses 
with durations between 1.3 and 200 ms. The range of 
pulse durations the analysis is sensitive to is reasonable 
given th e pulses observed from r otating radio transients 
(see, e.g. iMcLaughlin et al.ll2006l ). For each observation, 
single pulse events with signal-to- noise ratio (S/N) > 8 
were plotted and examined by eye. 

4. RESULTS 

None of the observations listed in Table [1] resulted in 
the detection of periodic or impulsive radio emission. 
The observations were used to place upper limits on any 
such emission, as we describe next. 

4.1. Upper Limits on Periodic Emission 

In order to compute upper limits on periodic emission 
for the observations analyzed in this work, we used the 
modified radiometer equation (e.g. iLorimer fc Kramed 
i20M), 



/3(S/N), 



[{Trcv + Tsky ) /G + Ssnr] 
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where Smin is the minimum detectable flux density in 
mjy, ji is the signal degradation factor due to quanti- 
zation, (S/N)^,-,j is the minimum signal-to-noise ratio 
considered, T^cv and T^^y are the receiver and sky tem- 
peratures in K, respectively, G is the telescope gain in K 
Jy~^, Ssnr is the flux density of the supernova remnant 
in Jy (if there is such an association), rip is the number 
of polarisations summed, tint is the integration time in 
seconds. A/ is the observing bandwidth in MHz, and S is 
the assumed duty cycle, ranging between and 1. Both 
the integration time and observing bandwidth were ad- 
justed downward to take into consideration data masked 
due to RFI removal. 

The duty cycle is related to the width according to 
S = Wb/P, however this width is not the intrinsic width 
of the pulsar's integrated profile. The intrinsic width is 
effectively broadened by the finite sampling time, t^amp, 
dispersive smearing within each channel, t£,M^ and multi- 
path scattering with the ISM, tscatt- The broadened 
width is related to the intrinsic width, Wi, according 
to 



i2 

^samp 



(2) 
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Here, the scattering time, tscatt, depends on the degree 
of inhomogeneity of the free electrons along the line-of- 
sight, which is also pre dicted using the NE2001 model 
(details can be found in lCordes fc La zio 2002). 

For the purposes of this work, the signal degradation 
factor due to quantization, /3, is ~ 1 since 16 bits are 
used. The gain of the GET, at 1950 MHz, ifO G = 
1.9 K Jy""'^. Also, the GBT's S-band receiver system is 
maintained at a temperature of Trcv = 20 K. The sky 
temperature, Tsky, is the sum of contributions due to 
Galactic sync h rotron radiation, as determined from the 
IHaslam et all (|1982[ ) all-sky radio mapQ, phis the 2.73 
K contribution from the cosmic microwave background. 
The supernova remnant associations and their respective 
flux densities are shown in Table [3] Finally, intrinsic 
pulse widths are taken to be between 3% and 50% of the 
pulse period. These widths are then broadened according 
to Equation O 

To determine the minimum detectable signal-to-noise 
ratio, (S/N)j^-„, in the presence of RFI, an observation of 
a known pulsar, PSR J1907-f0918, was used. Subsets of 
the observation of various durations were folded. The re- 
sulting folded profiles were examined by eye to determine 
an approximate threshold of detectibility, (S/N)^^^^^ = 4. 
Assuming purely white noise, a 4 cr detection has a prob- 
ability of P (4 cr) ~ 6.3 X 10~^ of occuring by chance. 

For magnetars for which a current ephemeris was avail- 
able (IE 1841-045, IE 2259+586 and 4U 0142+61; see 
Tabled, (S/N),^-^ — 4 was used since the exact period 
at the time of the observation was known. 

For sources for which there was an uncertainty on 
the folding frequency (SGR 1806-20, SGR 1900+14, 
and AX J1845-0258), the value of (S/N)^^„ used was 
larger, because of the larger number of Fourier bins 
searched. The number of bins searched was computed 
by estimating the frequency at the epoch of the obser- 
vation, which was extrapolated using a previously pub- 
lished frequency and frequency derivative, or in the case 
of AX J1845— 0258, a previously published frequency and 
an assumed magnetic field strength equal to 5 x 10^^ G 
(i.e. twice the largest inferred magnetic field strength 
for any known magnetar). A conservative fractional 
uncertainty on the frequency of lO"** was used for all 
three sources. This was done to take into considera- 
tion the possibility of a very large glitch, or anti-glitch. 
In all cases, the change in frequency due to a putative 
glitch was comparable to, or dominated by, the uncer- 
tainty in the extrapolated frequency. The frequency 
range searched was conservatively taken to be 3 times 
the change in frequency due to a potential glitch or anti- 
glitch, as described above. This range was divided by 
the frequency resolution, l/ti„t, to find the number of 
Fourier bins searched. This is equivalent to the num- 
ber of independent trials. In the case of SGR 1806—20 
and SGR 1900+14, the frequency range searched still 
amounted to only one Fourier bin. In the case of 
AX J1845-0258, 100 Fourier bins were searched. 

See the GBT prop oser's guide, 

|http://w ww.gb.nrao.edu/gbtprops/man/GBTpg.pdf' 

The electronic HEALPix version of the IHaslam et al.l 
iITQSI) 408 MHz map provided by NASA LAMBDA was used: 
http: / /lambda. gsfc.nasa.gov/product/forground/fg_haslam_get.cfm. 
The Galactic synchrotron e mission is assumed to have a power-law 
spectrum with index —2.5 IIKogut et al.ll20Tll 'l. 



For magnetar candidate GRB 050925, there is no in- 
formation about the spin period of the potential neu- 
tron star. Therefore, all Fourier bins were searched (i.e. 
5493164 and 4211426 bins for the observations on MJDs 
54053 and 54056, respectively). 

The probability corresponding to 4 cr was divided by 
the number of Fourier bins searched and converted back 
to an equivalent Gaussian sigma, which was used in 
Equation [TJ The values of (S/N)^-^^ used in Equation 
[T]are 5 for the observation of AX J1845-0258, and 6.79 
and 6.75 for the observations of GRB 050925 on MJDs 
54053 and 54056, respectively. 

For each observation, an upper limit on the fiirx den- 
sity of periodic emission has been computed as a function 
of pulse width. See Figure [T] for an example of the de- 
pendence of the upper limit on pulse width. Assuming a 
duty cycle of S = 5%, all observations have a minimum 
detectable fiux density of Smin ~ 0.013 mJy or lower. By 
assuming a distance for each source, this can be trans- 
lated to a minimum detectable luminosity of hmin ~ 1-6 
mJy kpc^ for all observations analysed. Table |4] summa- 
rizes the result^. 

4.2. Upper Limits on Single Pulses 

In addition to searching for periodic signals, the data 
were also searched for dispersed bright single pulses. 
Many single pulse events with signal-to-noise ratio > 8 
were detected, however none was found to be consistent 
with an astrophysical origin. 

To compute the minimum detectable flux density of 
single puls es in the observation s , the formalism intro- 
duced in Cordes & McLaughlin' (2003") was used. The 
relation between the brightness of a single pulse of astro- 
physical origin and its measured signal-to-noise ratio is 
given by 

S.J1I^.^. ,3, 

Wi y npAf 

where Si is the intrinsic flux density of a pulse of width 
Wi, (S/N)f, is the broadened signal-to-noise ratio, as mea- 
sured by matched filtering, Ssys — (Trcv + Tsky)/G + 
SsNB. is the system equivalent flux density, Wb is the 
broadened width of the single pulse. Wi and Wb are re- 
lated by Equation [21 In this work, the minimum signal- 
to-noise ratio considered for a single pulse is conserva- 
tively chosen to be (S/N)j,^-jj = 8. The large value of 
(S/N)j^-„ used is meant to exclude not only noise, but 
some RFI as well. For each observation, the minimum 
detectable luminosity of a 10-ms single pulse is reported 
in Table [H Examples of the dependence of Smin on pulse 
duration are shown in Figure [21 

Based on the previous discussion and the fact that no 
astrophysical single pulses were detected, we could con- 
clude that the sources observed do not emit single pulses 
brighter than Lmin. single, the minimum detectable sin- 
gle pulse luminosity. However, it is also possible the 
sources in question do occasionally emit single pulses 
bright enough to be detected, but the rates at which 

Note that the number of entries in Table[4]exceeds the number 
of entries listed in Table [l] because two observations were split due 
to RFI-induced jumps in the baseline. 
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these pulses are emitted are sufficiently small that no 
pulses were detected in our observations. With this in 
mind, a 3cr upper limit can be placed on the single pulse 
rate, knowing that < 1 single pulses were detected and 
assuming such pulses are emitted randomly. We have 
used the sum of all unmasked integration time for each 
source to compute limits on the single pulse rate. These 
apply only to pulses bright enough to be detected with 
durations between 0.33 ms and 200 ms, and are shown 
in Table El 

4.3. ToO Observations of 4U 0142+61 

ToO observations were triggered when an X-ray burst 
from 4 U 0142+61 was de t ected on MJD 54138 using 
RXTE (jGavriil et al.l [20071 [20TTI ). The first of six GBT 
ToO observations were on the day the burst was detected. 
Additional observations were scheduled at increasing in- 
tervals following the burst. The spacing between sub- 
sequent observations were: 1 day, 1 week, 2 weeks, 1 
month, and 6 months. None of these observations re- 
sulted in a detection of emission from 4U 0142+61. The 
collection of radio observations of 4U 0142+61, and a 
timeline of its activity in 2006 are shown in Figure [3] 

4.4. Constraints on Other Pulsars in the Field- of- View 

The surface density of known pulsars on the sky is large 
enough that a chance inclusion of an unrelated source in 
an observation is possible. We therefore also searched for 
signals at periods far from those of the magnetar targets. 

No new pulsars were discovered in any of the observa- 
tions searched. However, one previously known pulsar, 
PSR J19 07+0918, a 22 6-ms pulsar with a DM = 357 
cm~^ pc (jLorimer fc Xilo uris 2000), was detected in our 
only observation of SGR 1900+14, on MJD 54053. 

The minimum detectable flux density for an unknown 
pulsar as a function of spin period can be computed for 
the blind searches performed. Representative values are 
computed using Equationfl] assuming a duty cycle of (5 = 
5%, a period of P = 100 ms and DM = 100 cm'^ pc. 
The limits of the blind searches vary among targets since 
Tsky and Ssnr affect the system temperature. Also, the 
integration time and the RFI conditions (i.e. amount of 
data masked) are different for each observation. Typical 
results are reported in Table |4] and the dependences on 
DM and pulse period are shown in Figure S] 

5. DISCUSSION 

This work was designed to detect radio emission from a 
collection of magnetars and magnetar candidates at 1950 
MHz. The results presented here are more constraining 
than previous ly published limits at 1400 MHz for the 
same sources (|Burgav et aI1[2006L [Crawford et al.l [20071 : 
IGaensler et all [2005 Iden Hartog eFall [2007lf^ . How- 
ever, considering the variable nature of the three known 
radio-loud magnetars, the results of this work are also 
complimentary to the previously published constraints. 
Even though the variable nature of the radio-loud mag- 
netars limits the conclusions that can be drawn from in- 
dividual non-detections, it is still worthwhile to consider 
the ensemble of non-detections. 

^® Because of the large difference in observing frequencies, we 
do not directly compare our results with the limits obtained by 
ILorimer fc Xilouris! l|2000) for SGR 1900+14 at 430 MHz. 



The three known radio-loud magnetars have properties 
different from those of most other radio pulsars. Most 
importantly, the known radio magnetars are transient, 
i.e. their emission is often absent. Examples of this be- 
havior include the sudden appearance then later fading 
of XTE J181 0— 197's radio emiss ion over the course of 
~ 250 days ([Camilo et al.l l2007al ). as well as the more 
sporadic behavior observed in IE 1547.0—5408, which, 
for example, was not detected in observations on 2009 
January 22 and 2009 January 23, but was detected 
two days later on 2009 January 25 ( Camilo et al. 200!| 
IBurgav et al.l [20091 ). Also, when the known radio mag- 
netars are observed to be "on" , t heir luminosities are 
variable by a fac tor of a few (e.g. iCamilo et al.|[2007al : 
iLevin et al.ll2010i) . Therefore, the results presented in 
this work cannot be used to provide constraints when 
the sources were not being observed. 

5.1. Beaming Fractions of Magnetars 

Slowly rotating pulsars have relatively large light cylin- 
der radii and small polar caps, a nd thus are expected 
to ha ve narrow radio beams (e.g. iTauris fc Manchester 
119981 and references therein). Therefore, a magnetar 
may not be detected at radio frequencies because its ra- 
dio beam does not intersect the Earth. Unfortunately, 
it is difficult to estimate the fraction of the celestial 
sphere illuminated by any given pulsar without know- 
ing the geometry involved. The X-ray pulse profiles 
of m agnetars are typically broad ([Woods fc ThompsonI 
|2006[ ). thus their X-ray beaming fractions are likely quite 
large. Therefore, it is unreasonable to make assump- 
tions concerning the geometries of known magnetars 
based on their X-ray detections. Instead, naively as- 
suming that the radio beaming properties of magnetars 
are similar to those of rotation-powered radio pulsars, 
we use empirical results for the beaming fraction of ra- 
dio pulsars to attempt to constrain that of magnetars 
(Tauris & Manchester 1998). 

ITauris fc Manchesten (1998) estimate the fraction of 
pulsars, /, beamed towards the Earth as a function of 
spin period, P, to be 



/(P) = 0.09 



log 



P 



1 



0.03. 



(4) 



Given that there are three known radio-loud magne- 
tars out of a tota l of 13 sources observed (this work; 
Burgav et al.ll2006l: IGrawford et al.ll2007l: ILorimer et all 
2009f ) it is possible to estimate the probability that at 
least one of the undetected sources is beamed towards 
the Earth. Using a representative period of 7 s. Equa- 
tion 131 can be used to compute the probability that 4 
or more magnetars are beamed towards the Earth (that 
is, at least one of the undetected magnetars is beamed 
towards us). We find P(> 4) = 0.06%. This is equiva- 
lent to saying that if we had detected one of our targets, 
Equation|4|would be inconsistent with the magnetar data 
at the 3.2 cr level. Thus the non-detections reported here 
are consistent with being due to unfortunate beaming, 
assuming Gaussian statistics and that Equation [4| ap- 
plies. 

However, if we assume Equation [31 applies, and again 
consider representative period of 7 s, we find the proba- 
bility of three or more magnetars being beamed towards 
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us is P{> 3) ~ 0.75% (~ 2.4cr), which is suggestive 
that radio-magnetars have wider beams than rotation- 
powered pulsars with similar periods. In fact, by re- 
casting Equation |4] in terms of the pulsar beam's half- 
opening angle, p, we find that p ~ 30° is required 
for three out of 13 sources to be beamed towards us 
to be consistent with random orientations at the 3(t 
level. This value of p is ^ 10 times la rger than what 
is used by iTauris fc Manchesteii (|1998[ ). p - 2°, for 
P — 7 s, and is much more accommodating of the 
large duty cycles of the radio profiles of IE 1547.0—5408 
and PSR J1622-4950 (14% and ~1 1%, respectively; 
ICamilo et aLll2007bt Eevin et al.ll2010[ ). This suggested 
difference in p for magnetically powered and rotationally 
powered pulsars further reinforces the likely different ori- 
gin of the radio emission in these objects. 

6. CONCLUSIONS 

Using the GBT's Pulsar Spigot, five magnetars and two 
magnetar candidates were observed at 1950 MHz over the 
course of ^ 1 year. The data were searched for pulsed 
periodic emission, as well as bright single pulses. None of 
the targets observed were detected. The non-detections 
were used to place stringent constraints on the presence 
of any radio emission at the observing epochs. None 
of the seven targets observed showed periodic emission 
with luminosity larger than 1.6 mJy kpc^, or single pulse 
emission with luminosity larger than 7.6 Jy kpc^ at 1950 
MHz in any of the data sets. 

It is difficult to use the upper limits derived from this 
work to constrain the physical properties of the sources, 
since previous results suggest radio emission from mag- 
netars is highly variable. It is possible that observations 
were merely scheduled, by chance, at times when the 
magnetars' radio emission was too faint to be detected. 
On the other hand, it is also possible that the non- 
detections are nothing more than unfortunate beaming 
geometries. However, this work and previous work have 
combined to detect three out of 13 magnetars observed 
at radio frequencies. This suggests that the radio beam- 
ing fraction of magnetars is larger than that of rotation- 
powered pulsars with similar periods, in agreement with 
lines of reasoning. 

Finally, if any of the sources studied here becomes vis- 
ible at radio frequencies in the future, our results will be 
invaluable for providing a baseline with which detections 
can be compared. This may help us understand what 
magnetospheric conditions are required to produce radio 
emission, and the details of plasmas and currents in mag- 
netar magnetospheres. What is learned from magnetars 
could also likely help elucidate the radio emission mech- 
anism for pulsars in general. For this reason, the known 
magnetars should continue to be observed at radio fre- 
quencies. 
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Table 1 

Summary of GBT Observations of Five Magnetars and Two Magnetar Candidates. 



— n 

bourcc 


Uatc 


Epoch 


TTH 


ri„„i 

Dcci. 


Obs. length 






(MJD) 


(J2000) 


(J2000J 


(s) 
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o4uoz.yo 


1 fi- /I 1 ■ 1 o 


— U4.0D. 11 
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00. OZ.^ / 
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November 20, 2006 


54059.23 






10800 




iNOverriDer i^, zuud 


^AC\^'\ 1 H 
O^kjOo. io 




Dl ."rtO.UD 


9'^AO 

Z04:U 




November 17, 2006 


54056.71 






1500 




JNovember 20, zOOo 


54059.00 






16140 




reoruciry zuu/ 








o / ZU 




February 8, 2007 


541 -^Q S2^ 






3720 




February 15, 2007 


54146.33'' 






3600 




March 2, 2007 


54161. SI'' 






2400 




April 14, 2007 


54204.59=^ 






3180 




October 7, 2007 


54380.10'' 






6000 


AX J 1845-0258 


November 14, 2006 


54053.01 


18:44:55 


-02:56:56 


3600 




November 14, 2006 


54053.09 


20:13:47 


34:19:55 


3600 




November 17, 2006 


54056.73 






2760 


SGR 1806-20 


March 31, 2007 


54190.47 


18:08:40 


-20:24:37 


4500 




April 1, 2007 


54191.48 






3720 




April 21, 2007 


54211.47 






4260 


SGR 1900+14 


November 14, 2006 


54053.05 


19:07:15 


09:19:03 


3600 



Observation is part of Target-of-Opportunity obs ervation tr iggered after an X-ray burst 
was detected on MJD 54138 using RXTE tGavriil et al.ll20TllV 



Table 2 

Rotational Ephemerides Used for Folding. 



Parameter 


IE 1841-045 


IE 2259+586 


4U 0142+61 


l''requency (Hz) 

Frequency derivative (10^^^ Hz s^^) 

Epoch (MJD) 

Reference 


0.084863091(5) 
-289(2) 
54053 
Dib et al. (2008) 


0.14328613(6) 
-7(5) 
54050 
Dib (2009) 


0.1150918267(8) 
-26.89(6) 
53919 
Dib et al. (2007) 





Note. — Uncertainties reported are the 1 a uncertainties produced by TEMPO. Ephemerides 
are updated versions of what is reported in the references listed, and are valid only for short 
time intervals surrounding the observation epochs presented in this work. 



Table 3 

Summary of Relevant Properties of the Sources Observed. 



Source 


Dist.'' 


DM Max DM SNR Assoc. 


SSNR, 1950 


nn b 

sky 






Estimate Searched 






(kpc) 


(cm^'' pc) (cm"'' pc) 


(Jy) 


(K) 


Confirmed Magnetars 


IE 1841-045 


8.5 


800 2517 Kes 73 


3.8 


8 


IE 2259+586 


4.0 


150 1007 CTB 109 


15.8 


4 


4U 0142+61 


3.6 


100 503 




3.5 


SGR 1806-20 


8.7 


750 2517 




10 


SGR 1900+14 


15 


700 2014 




5.5 


Magnetar Candidates 


AX J 1845-0258 


8 


750 2517 G29. 6+0.1 


1.1 


8.5 


GRB 050925 


8.8 


300 1008 




4 



The distances reported here are the values used for estimating the DM, as well as for 
estimating luminosity limits later in this work. 

^ Sky temperatures include a 2.73 K contribution from the cosmic microwave backgrou nd, 
as wel l as a contribution from Galactic synchrotron radiation taken from the Haslam et al.l 
Ijl982l) all-sky 408 MHz map. Temperatures here are reported for 1950 MHz assuming a 
power-law spectrum with a synchrotron index of —2.5 IIKogut et al.|[20Tll ). 
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Table 4 

Upper Limits on Radio Emission from Magnetars and Magnetar Candidates. 



—5 

oource 
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r a b 
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— r '■ b c 
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(MJD) 


/ T \ 
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^ T 1 2 \ 
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(mJy) 


(Jy kpc^) 


(mJy) 


IE 1841—045 


54052.97 


0.0102 


0.74 


44.1 


3.2 


0.0173 


IE 2259+586 


54053.14 


0.0130 


0.21 
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1.1 
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0.0123 




04UD0. / i 


U.UilO 


U.lo 


Qn /I 
0U.4 
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0.0045 
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0.0078 
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0.0058 


0.08 
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0.4 


0.0099 
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n noQS 
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0.08 


32.9 


0.4 


0.0106 




54146.33 


0.0076 


0.10 


30.9 


0.4 


0.0128 




54161.51 


0.0078 


0.10 


32.8 


0.4 


0.0131 




54204.58 


0.0069 


0.09 


33.7 


0.4 


0.0117 




54380.10 


0.0048 


0.06 


32.5 


0.4 


0.0082 


AX J 1845-0258 


54053.01 


0.0092 


0.66 


38.4 


2.8 


0.0125 


GHB 050925 


54053.09 


0.0097 


0.75 


30.9 


2.4 


0.0097 




54056.73 


0.0130 


1.01 


33.3 


2.6 


0.0130 


SGR 1806-20 


54190.47 


0.0071 


0.54 


37.1 


2.8 


0.0121 




54191.48 


0.0078 


0.59 


37.2 


2.8 


0.0133 




54211.48 


0.0069 


0.52 


37.0 


2.8 


0.0117 


SGR 1900+14 


54053.05 


0.0071 


1.60 


33.8 


7.6 


0.0121 



^ A 5% duty cycle is assumed. 

Distances assumed are listed in Tabled 

The values reported here are for a duration of 10 ms. Single pulse durations searched are in the 
range 0.33 ms to 200 ms. 

Blind periodicity search. Assumed 5% duty cycle, 100 ms period, and DM = 100 cm~^ pc. 

The observation was split into smaller portions that were searched independently due to significant 
changes in the data mean caused by RFI (see §3). 



Table 5 

Upper Limits on Single Pulse Rates of Magnetars and Magnetar Candidates. 



Source 


Total usable time 


Minimum detectable 


Single pulse rate limit 




(hr) 


single pulse luminosity, -L1950 
(Jy kpc2) 


(hr-l) 


IE 1841-045 


0.66 


3.2 


4.5 


IE 2259+586 


3.5 


1.1 


0.85 


4U 0142+61 


8.5 


0.4 


0.35 


AX J1845-0258 


0.98 


2.8 


3.1 


GRB 050925 


1.6 


2.4 


1.9 


SGR 1806-20 


2.8 


2.8 


1.1 


SGR 1900+14 


0.74 


7.6 


4.1 



Note. — The single pulse rate limit applies to pulses brighter than the limiting L1950 and 
with durations between 0.33 ms and 200 ms. 



10 




0.1 0.2 0.3 0.4 0.5 

Duty Cycle (%) 



Figure 1. Best flux density limits at 1950 MHz as a function of duty cycle for periodic emission from each of the seven targets observed. 
From top to bottom the curves correspond to IE 2259+586, IE 1841-045, GRB 050925 (dashed), AX J1845-0258, SGR 1900+14, 
SCR 1806-20, and 4U 0142+61. 




10° 10^ 10^ 

Pulse width (bins) 



Figure 2. The best flux density limits for single pulses at 1950 MHz as a function of pulse width for each of the seven targets observed. 
From top to bottom the curves correspond to IE 2259+586, IE 1841-045, AX J1845-0258, SGR 1806-20, SGR 1900+14, GRB 050925 
(dashed), and 4U 0142+61. 
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Figure 3. A timeline of activity f or 4U 0142+61 including a glitch (dotted line), six X-ray bursts at three different epochs (dashed 
lines, from lGavriil et al.|[2007l. 120111) . and upper limits on radio emission at 1950 M Hz (downward-poin t ing triangles). Note that typical 
luminosities for the radio-detected magnetars are significantly off the top of the plot IjCamilo et al.ll2006l . I2007bl : ILevin et al.ll2010l ). 
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Figure 4. Typical flux density limits for a blind search for pulsars as a function of spin period for various DMs. These limits are for the 
observation of AX J1845— 0258 on MJD 54053. An intrinsic duty cycle of 5% is assumed. The duty cycle is broadened by the observing 
system's finite sample time, dispersive smearing and scattering. From left to right the lines correspond to dispersion measures of cm~^ 
pc, 500 cm~'^ pc, 1000 cm~^ pc, and 2000 cm~^ pc. 



